Ab initio calculations of the atomic and electronic structure of CaTiOs (001) and (Oil) 

surfaces 
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We present the results of calculations of surface relaxations, energetics, and bonding properties 
for CaTiOa (001) and (Oil) surfaces using a hybrid B3PW description of exchange and correlation. 
We consider both CaO and Ti02 terminations of the non-polar (001) surface, and Ca, TiO and O 
terminations of the polar (Oil) surface. On the (001) surfaces, we find that all upper-layer atoms 
relajc inwards on the CaO-terminated surface, while outward relaxations of all atoms in the second 
layer are found for both terminations. For the Ti02-terminated (001) surface, the largest relaxations 
are on the second-layer atoms. The surface rumpling is much larger for the CaO-terminated than 
for the Ti02-terminated (001) surface, but their surface energies are quite similar at 0.94 eV and 
1.13 eV respectively. In contrast, different terminations of the (Oil) CaTiOs surface lead to very 
different surface energies of 1.86 eV, 1.91 eV, and 3.13 eV for the 0-terminated, Ca-terminated, and 
TiO-terminated (Oil) surface respectively. Our results for surface energies contrast sharply with 
those of Zhang et al. [Phys. Rev. B 76, 115426 (2007)], where the authors found a rather different 
pattern of surface energies. We predict a considerable increase of the Ti-O chemical bond covalency 
near the (Oil) surface as compared both to the bulk and to the (001) surface. 

PACS numbers: 68.35.Ct, 68.35.Md, 68.47.Gh 



o 



0^ 

O 

00 

O 



X 



I. INTRODUCTION 

Oxide perovskites are in demand for a variety of in- 
dustrial applications as a result of their diverse physi- 
cal properties i^i^ For example, CaTiOs is a cubic per- 
ovskite that is widely used in electronic ceramic ma- 
terials and as a key component of synthetic rock to 
immobilize high-level radioactive waste.^ Thin films of 
ABO3 perovskite ferroelectrics are important for many 
applications In particular, the titanates are interesting 
materials regarding their electrochemical properties and 
are promising as components for electrodes and sensors. 
Surface properties of CaTiOs are important for cataly- 
sis and for epitaxial growth of high Tc superconductors. 
For all these applications, the surface structure and the 
associated surface electronic and chemical properties are 
of primary importance. 

In view of this technological importance, it is surprising 
that there have been so few ah initio studies of CaTiOs 
surface atomic and electronic structure. For the CaTiOs 
(001) surface we are only aware of the work of Wang et 
al^ and Zhang et al.^ In contrast, several other ABO3 
perovskite (001) surfaces have been widely studied. For 
example, ah and 
classical shell-mode P^'^'^ studies were published for the 
(001) surfaces of SrTiOs. The (GDI) surfaces of cu- 
bic perovskites have also been extensively investigated 
experimentally. For example, the SrTiOs (001) sur- 
face relaxations and rumplings have been studied by 
means of low energy electron diffraction (LEED), reflec- 
tion high-energy electron diffraction (RHEED), medium 
energy ion scattering (MEIS), and surface x-ray diffrac- 
tion (SXRD) measurements The status of the 
degree of agreement between theory and experiment for 



these SrTiOa surfaces is summarized in Ref. Q- 

ABO3 perovskite (Oil) surfaces are considerably less- 
well studied than (001) surfaces, both experimentally 
and theoretically. However, there has been a surge of 
recent interest, focusing mainly on SrTiOs, in which 
STM, UPS, XPSr^i^i and Auger spectroscopies as well 
as LEED studies^^i^^i^i^^i^ have been carried out. On 
the theory side, the first ah initio calculations for SrTiOs 
(Oil) surfaces were performed by Bottin et al.,^^ who 
carried out a systematic study of the electronic and 
atomic structures of several (1x1) terminations of the 
(Oil) polar orientation of the SrTiOs surface. They 
found that the electronic structure of the stoichiomet- 
ric SrTiO and O2 terminations showed marked differ- 
ences with respect to the bulk as a consequence of the 
polarity compensation. Later, Heifets et al.^^ performed 
ah initio Hartree-Fock (HE) calculations for four possi- 
ble non-polar terminations (TiO, Sr, and two kinds O 
terminations) of the SrTiOs (Oil) surface. The authors 
found that the surface energy of the O-terminated (Oil) 
surface is close to that of the (001) surface, suggesting 
that both (Oil) and (GDI) surfaces can coexist in poly- 
crystalline SrTiOa. Most recently, we performed an ab 
initio study of SrTiOa (Oil) surfaces^ using a hybrid 
Hartree-Fock (HE) and density-functional theory (DET) 
exchange-correlation functional, in which HE exchange 
is mixed with Becke's three-parameter DET exchange 
and combined with the nonlocal correlation functional of 
Perdew and Wang (B3PW)42ii^ Our calculations indi- 
cated a remarkably large increase in the Ti-0 bond cova- 
lency at the TiO-terminated (Oil) surface, significantly 
larger than for the (GGl) surfaces. 

Regarding other ABO3 (Oil) surfaces, Heifets et al^ 
investigated the atomic structure and charge redistribu- 
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tion of different terminations of BaZrOs (Oil) surfaces 
using density- functional methods. They found that the 
0-terminated (Oil) surface had the smallest cleavage en- 
ergy among (Oil) surfaces, but that this value was still 
twice as large as the cleavage energy needed for the for- 
mation of a pair of complementary (001) surfaces. More- 
over, we recently performed ah initio B3PW calculations 
for the technologically important BaTiOa and PbTiOs 
(Oil) surfacesi^ Our calculated surface energies showed 
that the Ti02-terminated (001) surface is slightly more 
stable than the BaO- or PbO-terminated (001) surface 
for both materials, and that 0-terminated BaTiOs and 
TiO-terminated PbTiOa (Oil) surfaces have surface en- 
ergies close to that of the (001) surface. 

The only existing ah initio study of CaTiOa (Oil) polar 
surfaces was performed by Zhang et alM^ In addition to 
the (001) surfaces, they studied four possible non-polar 
terminations of the (Oil) surface, namely the TiO, Ca, 
asymmetric A- type O, and symmetric B-type O termina- 
tions. The results indicated that the most favorable sur- 
faces are the CaO-terminated (001) surface, the A-type 
0-terminated (Oil) surface, and the Ti02-terminated 
(001) surface, in that order. 

With the sole exception of the calculation on 
CaTiOs by Zhang et al.^ all of the first- principles 
and shell-model studies of ABO3 perovskite surface 
energie o'''-^^i-^^'^'^i'^°'^^i'^^ have found that the lowest- 
energy (001) surface is lower in energy than any of the 
(Oil) terminations. Zhang et al.^ on the contrary, re- 
ported a surface energy of 0.837 eV for their "A-type" 
0-terminated (Oil) surface of CaTiOs, to be compared 
with 1.021 eV for the Ti02-terminated (001) surface. Be- 
cause this result contrasts sharply with the other previous 
calculations, we were particularly motivated to check this 
result independently in our current study. 

In this study, we have performed predictive ah initio 
calculations for CaTiOa (001) and (Oil) surfaces, using 
the same B3PW approach as in our previous work.^''^'^ 
As in the work of Zhang et al.^ we do not explicitly 
include octahedral rotations in the surface calculations, 
even though such rotations are likely to be more impor- 
tant for CaTiOs than for many other perovskites; we dis- 
cuss and justify this approximation at the end of Sec. II. 
In contradiction to the work of Zhang et al.^ we find 
that the pattern of surface energies of CaTiOa is similar 
to that of other perovskites. In particular, we find that 
the 0-terminated CaTiOs (Oil) surface is higher in en- 
ergy than either of the Ti02- or CaO-terminated (001) 
surfaces. We also report the surface relaxations and rum- 
plings and the charge redistributions and changes in bond 
strength that occur at the surface. 

The manuscript is organized as follows. In Sec. II we 
present our computational method and provide details 
of the surface slab models on which the calculations were 
performed. The results of our calculations for surface 
structures, energies, charge distributions, and bond pop- 
ulations are reported in Sec. III. Finally, we discuss the 
results and present our conclusions in Sec. IV. 



II. COMPUTATIONAL METHOD AND 
SURFACE SLAB CONSTRUCTION 

To perform the first-principles DFT-B3PW calcula- 
tions we used the CRYSTAL-2003 computer code^i^ 
which employs Gaussian-type functions (GTFs) localized 
at atoms as the basis for an expansion of the crystalline 
orbitals. The features of the CRYSTAL-2003 code that 
are most important for this study are its ability to cal- 
culate the electronic structure of materials within both 
Hartree-Fock and Kohn-Sham Hamiltonians and its abil- 
ity to treat isolated 2D slabs without artificial repetition 
along the z-axis. However, in order to employ the linear 
combination of atomic orbitals (LCAO)-GTF method, it 
is desirable to have optimized basis sets (BS). The BS op- 
timization for SrTiOa, BaTiOa, and PbTiOa perovskites 
was developed and discussed in Ref. Here we em- 

ploy this BS, which differs from that used in Refs. [T^ITsj 
by inclusion of polarizable d-orbitals on O ions. It was 
shown^ that this leads to better agreement of the cal- 
culated lattice constant and bulk modulus with experi- 
mental data. For the Ca atom we used the same BS as 
in Ref. [45|. 

Our calculations were performed using the hybrid 
exchange-correlation B3PW functional involving a hy- 
brid of non-local Fock exact exchange, LDA exchange and 
Becke's gradient corrected exchange functional com- 
bined with the nonlocal gradient corrected correlation 
potential by Perdew and Wang."*^ The Hay-Wadt small- 
core effective core pseudopotentials (ECP) were adopted 
for Ca and Ti atoms4^ The small-core ECP's replace 
only the inner core orbitals, while orbitals for sub- valence 
electrons as well as for valence electrons are calculated 
self-consistently. Oxygen atoms were treated with the 
all-electron BS. 

The reciprocal space integration was performed by 
sampling the Brillouin zone of the five-atom cubic unit 
cell with an 8 X 8 X 8 Pack-Monkhorst grid for the bulk, — 
and an 8x8 grid for the slab structure, providing a bal- 
anced summation in direct and reciprocal spaces. To 
achieve high accuracy, large enough tolerances of 7, 8, 7, 
7, and 14 were chosen for the Coulomb overlap. Coulomb 
penetration, exchange overlap, first exchange pseudo- 
overlap, and second exchange pseudo-overlap parame- 
ters, respectively!^ 

The CaTiOa (001) surfaces were modeled with two- 
dimensional slabs consisting of several planes perpendic- 
ular to the [001] crystal direction. To simulate CaTiOs 
(001) surfaces, we used slabs consisting of seven alter- 
nating Ti02 and CaO layers, with mirror symmetry pre- 
served relative to the central layer. The 17-atom slab 
with CaO-terminated surfaces and the 18-atom slab with 
Ti02-terminated surfaces are shown in Figs. 1(a) and (b) 
respectively. These slabs are non-stoichiometric, with 
unit-cell formulae Ca4Ti30io and Ca3Ti40ii, respec- 
tively. These two (CaO and Ti02) terminations are the 
only possible flat and dense (001) surface terminations 
of the perovskite structure. The sequence of layers with 
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FIG. 1: (Color online.) Side view of CaTiOs (001) surfaces, 
(a) CaO-terminated surface, (b) Ti02-terminated surface, 
with definitions of surface rumpling s and the near-surface 
interplanar separations Adi2 and Ad23. 



(001) orientation, and the definitions of the surface rum- 
pling s and the interplane distances Ac?i2 and Ac?23, are 
illustrated in Fig. 1. 

The problem in modeling the CaTiOs (Oil) polar sur- 
face is that, unlike the CaTiOs (001) neutral surface, it 
consists of charged 0-0 and CaTiO planes, as illustrated 
in Fig. 2. Assuming nominal ionic charges of 0^~, Ti*+, 
and Ca^+, a simple cleavage would create a negatively- 
charged 0-0 surface and a positively-charged CaTiO 
surface, leading either to an infinite macroscopic dipole 
moment perpendicular to the surface for a stoichiomet- 
ric slab terminated by planes of different kinds (O2 and 
CaTiO) as in Fig. 3(a), or to a net infinite charge for a 
non-stoichiometric symmetric slab as shown in Figs. 3(b) 
and (c). It is known that such crystal terminations make 
the surface unstable. In proper first-principles calcu- 
lations on slabs of finite thickness, charge redistributions 
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FIG. 2: (Color online.) Sketch of the cubic CaTiOs perovskite 
structure showing two (Oil) cleavage planes that give rise to 
charged CaTiO and O2 (Oil) surfaces. 




FIG. 3: (Color online.) Possible (Oil) surface slab mod- 
els considered in the text, [(a)-(c)] Slabs obtained by sim- 
ple cleavage, yielding mixed, O-terminated, and CaTiO- 
terminated polar surfaces, respectively, [(d)-(f)] Slabs with 
nonpolar TiO-terminated, Ca-terminated, and O-terminated 
surfaces, respectively. 



near the surface arising during the self-consistent field 
procedure could, in principle, compensate at least par- 
tially for these effects. However, previous careful studies 
for SrTiOg^i^i^ have demonstrated that the resulting 
surfaces have a high energy, and that the introduction of 
surface vacancies provides an energetically less expensive 
mechanism for compensating the surfaces. 

For these reasons, we limit ourselves here to non-polar 
CaTiOs (Oil) surfaces that have been constructed by 
modifying the composition of the surface layer. Remov- 
ing the Ca atom from the upper and lower layers of 
the 7-layer symmetric CaTiO-terminated slab generates 
a neutral and symmetric 16-atom supercell with TiO- 
terminated surfaces as illustrated in Fig. 3(d). Remov- 
ing both the Ti and O atoms from the upper and lower 
layers of the 7-layer symmetric CaTiO-terminated slab 
yields a neutral and symmetric 14-atom supercell with 
Ca-terminated surfaces as shown in Fig. 3(e). Finally, 
removing the O atom from the upper and lower layers 
of the 7-layer symmetric 0-0 terminated slab, we ob- 
tain the neutral and symmetric 15-atom supercell with 
O-terminated surfaces shown in Fig. 3(f). The stoichiom- 
etry of these surface terminations, and the number of 
bonds cleaved, are comprehensively discussed for the case 
of SrTiOs in Ref. 

Before leaving this Section, it is worth discussing the 
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TABLE I; Calculated effective charges Q and bond popula- 
tions P (in e) for bulk CaTiOa. 



Ion or bond Property Value 

Ca Q 1.782 

O Q -1.371 

Ti Q 2.330 

Ca-O P 0.006 

Ti-0 P 0.084 

O-O P -0.010 



TABLE II: Computed atomic relaxation (in percent of the 
bulk lattice constant ao) for the Ti02- and CaO-terminated 
CaTiOa (001) surfaces. Positive values indicate outward dis- 
placements. 



CaO-terminated Ti02-terminated 
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This study 
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This study 
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-1.71 
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Ti 


1.12 


Ca 


2.75 




O 


0.01 





1.05 



issue of the tilting of TiOg octahedra in CaTiOa. Xray 
and neutron diffraction studies have not definitively es- 
tablished the phase-transition sequence at higher tem- 
perature, but clearly show that the crystal adopts an 
orthorhombic structure with space group Pbnm below 
~1380K.^i^i^'Si This room-temperature ground state 
has a 20-atom unit cell and is a slight modification from 
the ideal perovskite structure involving a pattern of tilts 
of the TiOg octahedra according to the a~ a~ pattern 
in Glazer's notation.^ The octahedral tilts have also 
been studied using first-principles calculations Be- 
cause these tilts are substantial (^10°),^'' it is possible 
that they may have some impact on the surface structure 
and energetics. However, we have not included octahe- 
dral tilts in the work presented here for several reasons. 
First, we want to compare with previous calculations, 
which have universally not included octahedral tilts. Sec- 
ond, the CRYSTAL-2003 code package does not provide 
for efficient structural optimization as would be needed 
to study these tilts, and the larger surface unit cells that 
would be required would make the calculations impracti- 
cal. But finally and most importantly, we estimate that 
the energy scale of the tilts (^0.1 eV) is small compared 
to the energy scale of the surface cleavage and relaxation 
energies (a few eV), so that it is reasonable to neglect 
them in a first study. The interactions between bulk tilts 
and surface relaxations remains an interesting question 
for future study. 

III. RESULTS OF CALCULATIONS 

A. CaTiOs bulk atomic and electronic structure 

As a starting point of our calculations, we calculated 
the CaTiOa bulk lattice constant and found it to be 
3.851 A, slightly smaller than the experimental value of 
3.895 Ai ^^i^^ We used the theoretical bulk lattice constant 
in the following surface structure calculations. To char- 
acterize the chemical bonding and covalency effects, we 
used a standard MuUiken population analysis for the ef- 
fective atomic charges Q, bond populations P, and other 
local properties of the electronic structure as described 
in, e.g., Refs. fsollsoj . Our calculated effective charges 
and bond populations for bulk CaTiOs are presented in 



Table I. The bond population of the Ti-0 bond is clearly 
much larger than that of the Ca-O bond, consistent with 
partial Ti-0 covalency, and the small but negative 0-0 
population indicates a repulsive overlap of oxygen shells 
in bulk CaTiOs. 



B. CaTiOs (001) surface structure 

The atomic displacements obtained using the ah initio 
B3PW method for Ti02- and CaO-terminated CaTiOs 
(001) surfaces are shown in Table II. According to the re- 
sults of our calculations, atoms of the first surface layer 
relax inwards, i.e. towards the bulk, for both Ti02 and 
CaO-terminated (001) surfaces. The latter result is in 
disagreement with the previous calculations of Wang et 
al.,^ who calculated that the first-layer oxygen atoms on 
the (001) surface should relax outwards by 0.7% of the 
bulk lattice constant oq. According to our calculations, 
they move inwards by 0.42% of oq. Our calculated in- 
ward relaxation of the first-layer oxygen atoms on the 
CaO-terminated CaTiOa (001) surface is in line with pre- 
vious ab initio studies dealing with BaTiOa, PbTiOa, 
and BaZrOa (001) surfacesj^Siiai^iS but contrasts with 
the outward relaxation of first-layer oxygen atoms on the 
SrO-terminated SrTiOa (001) surface i^^i^ According to 
the results of our current calculations, outward relax- 
ations are found for all atoms in the second layer for 
both CaO and Ti02 terminations of the CaTiOa (001) 
surface. 

Table II shows that the relaxations of the surface metal 
atoms are much larger than those of the oxygens on both 
the Ti02- and CaO-terminated CaTiOa (001) surfaces, 
leading to a considerable rumpling of the outermost sur- 
face plane. For the Ti02-terminated case, we found much 
larger displacements in the second layer than in the first 
layer. This behavior contrasts with the atomic relaxation 
pattern of the Ti02-terminated BaTiOa (001) surface, 
where the upper-layer Ti relaxation is generally larger 
than the second-layer Ba relaxationi^2i^ However, it is 
in line with the only existing ab initio study of the Ti02- 
terminated CaTiOa (001) surface,^ as well as with other 
ab initio studies dealing with related B02-terminated 
ABOa (001) surfaces, such as PbTi0332-62 and BaZrOa^ 
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TABLE III: Calculated surface rumpling s, and relative dis- 
placements Adij between the three near-surface planes, for 
the CaO- and Ti02-terminated CaTiOs (001) surface. Units 
are percent of the bulk lattice constant. 

CaO-terminated Ti02-terminated 

S Adi2 Ad23 S Adi2 Ad23 

7.89 -9.43 1.12 1.61 -4.46 2.75 



where the second-layer anion (Pb or Ba) relaxations were 
larger than the upper-layer (Ti or Zr) ones. 

In order to compare the calculated surface structures 
with experimental results, the surface rumpling s (the 
relative oxygen displacement relative to the metal atom 
in the surface layer) and the changes in interlayer dis- 
tances Adi2 and Ac?23 (where 1, 2 and 3 label the near- 
surface layers) are presented in Table III. Our calcula- 
tions of the interlayer distances are based on the posi- 
tions of the relaxed metal ions (Fig. 1), which are known 
to be much stronger electron scatters than the oxygen 
ionsi^ The amplitude of the surface rumpling on the 
CaO-terminated surface is predicted to be almost five 
times larger than that for the Ti02-terminated (001) one. 
From Table III, one can see that both CaTiOa (001) sur- 
faces show a reduction of interlayer distance Adi2 and an 
expansion of Ad23. The reduction of interlayer distance 
Adi2 is twice as large for the CaO-terminated surface 
than it is for the Ti02-terminated surface. Our calcu- 
lations dealing with the surface rumpling s, reduction of 
interlayer distances Adi2, and expansion of interlayer dis- 
tances Ad23 are in qualitative agreement with the only 
existing ab initio study dealing with CaTiOa (001) sur- 
face structuresj^ 

To the best of our knowledge there are no experimental 
measurements with which we can compare our caculated 
values of s, Adu, and Ad23 on the CaTiOs (001) sur- 
faces. Even when such data do exist, it is sometimes con- 
tradictory, as is the case for the SrO-terminated SrTiOs 
(001) surface, where existing LEED^i and RHEED25 ex- 
periments contradict each other regarding the sign of 
Adi2. 

The calculated atomic displacements, effective static 
charges, and bond populations between nearest metal 
and oxygen atoms are given for the Ti02- and CaO- 
terminated (001) surfaces in Table IV. The major effect 
observed here is a strengthening of the Ti-0 chemical 
bond near the surface. Recall from Table I that the 
Ti and O effective charges (2.330 e and — 1.371 e, respec- 
tively) in bulk CaTiOs are much smaller than expected 
in an ideal ionic model, and that the Ti-0 bond pop- 
ulation is 0.084 e. Table IV shows that the Ti-0 bond 
population for the Ti02-terminated (001) surface is con- 
siderably larger than the associated bulk value. Compar- 
ing with the very small bulk Ca-0 bond populations of 
0.006 e from Table I, we see that the Ca-0 bond popu- 
lations near the CaO-terminated (001) surface in Table 
IV are more than three times larger than in the bulk. 



TABLE IV: Calculated absolute magnitudes of atomic dis- 
placements D (in A), effective atomic charges Q (in e), and 
bond populations P (in e) between nearest metal-oxygen 
pairs, for the for the Ti02- and CaO-terminated CaTiOs 
(001) surfaces. 

Layer Property Ion Ti02-terminated Ion CaO-terminated 
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but more than five times smaller than the Ti-0 bond 
populations on the Ti02-terminated (001) surface. 

C. CaTiOa (Oil) surface structures 

As explained in Sec. II, non-polar TiO-, Ca-, and O- 
terminated surfaces can be constructed for the CaTiOs 
(Oil) surface as in Figs. 3(d)-(f) respectively. Details of 
the relaxed structures obtained from our calculations for 
these three terminations are given in Tables V and VI. 

On the TiO-terminated (Oil) surface, the upper-layer 
Ti atoms move inwards by 7.14% of the bulk lattice 
constant ao, whereas the O atoms move outwards by 
4.67% (Table V), leading to a large surface rumphng of 
11.81% (Table VI), in excellent agreement with the cor- 
responding surface rumpling of 12.10% calculated earlier 
by Zhang et al^ The second-layer oxygen atoms move 
inwards by less than 1% of oq. The displacement mag- 
nitudes of the atoms in the third layer are larger than 
in the second layer, but smaller than in the top layer. 
The Adi2 values in Table VI show that the reduction of 
the distance between the first and second layers is three 
times larger than the corresponding expansion between 
the second and third layers. 

On the Ca-terminated (Oil) surface. Table V shows 
that the Ca atoms in the top layer move inwards very 
strongly, while the O atoms in the second layer only move 
outwards very weakly. The pattern of oxygen displace- 
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TABLE V; Calculated atomic relaxations of the CaTiOa (Oil) 
surfaces (in percent of the bulk lattice constant ao) for the 
three surface terminations. Positive signs correspond to out- 
ward displacements. 



Layer Ion Az Ay 

TiO-terminated surface 

1 Ti -7.14 

1 O 4.67 

2 O -0.44 

3 Ca -2.75 
3 O -3.79 
3 Ti -0.78 

Ca-terminated surface 

1 Ca -16.05 

2 O 1.35 

3 Ti -0.37 
3 O -1.71 
3 Ca -0.93 

O-terminated surface 

1 O -6.10 -2.16 

2 Ti -0.26 -4.70 
2 Ca -2.10 -0.27 

2 O 3.43 8.05 

3 O -0.55 1.90 



ments is similar to that found on the TiO-terminated 
(Oil) surface, in that the inward oxygen displacement in 
the third layer is larger than the outward displacement 
in the second layer, but the Ti and Ca displacements in 
the third layer are smaller than the second-layer oxygen 
atom displacements. 

The O-terminated (Oil) surface has sufficiently low 
symmetry that some displacements occur in the y as 
well as in the z direction. The O atoms in the top 
layer move mostly inwards (~6%) but also have some dis- 
placement along the surface (^2%). On the other hand, 
the second-layer Ti atoms on this surface move strongly 
along the surface, and also slightly inwards. The second- 
layer Ca atoms move slightly in the same y direction, and 
also inwards, while the second-layer O atoms move very 
strongly in the y direction (but in the opposite direction 
compared to the top-layer O atoms) and rather strongly 
outwards. The third-layer O atoms move in the same di- 
rection as the second-layer O atoms along the y-axis, but 
their displacement magnitude are more than four times 
smaller, and they also move slightly inwards. Table VI 
shows that there is a substantial contraction of the inter- 
layer distance Adi2 and only a very slight expansion of 



TABLE VL Surface rumpling s and relative displacements 
Adij (in percent of the bulk lattice constant ao) for the three 
near-surface planes on the TiO- and O-terminated CaTiOs 
(Oil) surfaces. 



TiO terminated 




O terminated 


s Adi2 


Ad23 


Adi2 Ad23 


11.81 -6.70 


2.31 


-5.84 0.29 



D. CaTiOa (001) and (Oil) surface energies 

In the present work, we define the unrelaxed surface 
energy of a given surface termination A to be one-half 
of the energy needed to cleave the crystal rigidly into an 
unrelaxed surface A and an unrelaxed surface with the 
complementary termination A'. For CaTiOa, for exam- 
ple, the unrelaxed surface energies of the complementary 
CaO- and Ti02-terminated (001) surfaces are equal, as 
are those of the TiO- and Ca-terminated (Oil) surfaces. 
The relaxed surface energy is defined to be the energy of 
the unrelaxed surface plus the (negative) surface relax- 
ation energy. These definitions are chosen for consist ency 
with Refs. [Blli. Unlike the authors of Refs. [UllllfiJ, 
we have made no effort to introduce chemical potentials 
here. Thus, while the values of the surface energies E'surf 
reflect the cleavage energies and thus give some informa- 
tion about trends in the surface energetics, they should 
be used with caution when addressing questions of the 
relative stability of surfaces with different stoichiometrics 
in specific environmental conditions. 

To calculate the CaTiOs (001) surface energies, we 
start with the cleavage energy for the unrelaxed CaO- 
and Ti02-terminated surfaces. In our calculations the 
two 7-layer CaO- and Ti02-terminated slabs, containing 
17 and 18 atoms respectively, represent together 7 bulk 
unit cells of 5 atoms each. Surfaces with both termina- 
tions arise simultaneously under cleavage. According to 
our definition, we assume that the relevant cleavage en- 
ergy is distributed equally between created surfaces, so 
that both the CaO- and Ti02-terminated surfaces end 
up with the same unrelaxed surface energy 

^LT/' = Jl^r '(CaO) + 4T^(Ti02) - 7i?,.,k], (1) 

where E!,^^^\Ga.O) and i?g|!^"^(Ti02) are the unrelaxed 
CaO- and Ti02-terminated slab energies, -Ebuik is the en- 
ergy per bulk unit cell, and the factor of 1/4 comes from 
the fact that we create four surfaces upon the cleavage 
procedure. Our calculated unrelaxed surface energy for 
these surfaces is 1.40 eV, as shown in Table VII. The cor- 
responding relaxation energies are calculated using 

i?(,,,)(A) = i[41'(A)-4x'(A)], (2) 

where A = CaO or Ti02 and E^l^^i-^) is the slab energy 
after both sides of the slab have been allowed to relax. 
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TABLE VII: Calculated cleavage, relaxation, and surface en- 



ergies for CaTiO 
cell). 


3 (001) and (Oil) surfaces (in 


eV per 


surface 


Surface 


Termination 


suri 


Erel 


Egurf 


CaTiOs (001) 


Ti02 


1.40 


-0.27 


1.13 




CaO 


1.40 


-0.46 


0.94 


CaTiOs (Oil) 


TiO 


4.61 


-1.48 


3.13 




Ca 


4.61 


-2.70 


1.91 




O 


:',.:!() 


-i.ll 


1.86 



where i?surf(0) and -£'j,j™J(0) arc the surface energy and 
the relaxed slab total energy for the O-terminated (Oil) 
surface. The results are again summarized in Table VII. 
Unlike for the (001) surface, we sec that different ter- 
minations of the (Oil) surface lead to large differences 
in the surface energies. Here the lowest calculated sur- 
face energy is 1.86 eV for the O-terminated (Oil) surface, 
while the TiO-terminated (3.13 eV) is much more costly 
than the Ca-terminated (Oil) surface (1.91 eV). 



We find relaxation energies of —0.27 and — 0.46eV for 
the Ti02-terminated and CaO-terminated surfaces, re- 
spectively. The final surface energies are then obtained 
as a sum of the cleavage and relaxation energies using 

£;,urf(A) = EiZ\^) + ^(rel)(A). (3) 

The resulting surface energies of the two (001) surfaces 
arc comparable, but that of the Ti02-terminatcd surface 
is slightly larger than that of the CaO-terminated one 
(1.13 vs. 0.94eV), as summarized in Table VII. 

In order to calculate the surface energies of the TiO- 
and Ca-terminated surfaces shown in Fig. 3(d) and (e), 
containing 16 and 14 atoms respectively, we start with 
the cleavage energy for unrelaxed surfaces. The two 7- 
plane Ca- and TiO-terminated slabs represent together 
six bulk unit cells. The surfaces with both terminations 
arise simultaneously under cleavage of the crystal, and 
the relevant cleavage energy is divided equally between 
these two surfaces, so we obtain cleavage energies accord- 
ing to 

E^:Pi^) = \[E^^\Ga) + 4:b^(TiO) - 6^;^.,^ (4) 

where A denotes Ca or TiO, £^sia"''(A) is the energy of 
the unrelax;ed Ca or TiO terminated (Oil) slab, £^buik 
is the energy per bulk unit cell, and again the factor of 
1/4 arises because four surfaces arc created upon cleav- 
age. Our calculated cleavage energy for the Ca or TiO- 
terminated (Oil) surfaces of 4.61 eV is more than three 
times larger than the relevant cleavage energy for the 
CaO- or Ti02-tcrminated (001) surfaces. Finally, the 
surface energy Es„^{{A) is just a sum of the cleavage and 
relaxation energies, as in Eq. (3). 

When wc cleave the crystal along (Oil) in another way, 
as in Fig. 3(f), wc obtain two identical O-terminated sur- 
face slabs containing 15 atoms. The cleavage energy of 
3.30 eV computed for this O-terminated surface is slightly 
smaller than for the Ca or TiO-terminated (Oil) surfaces, 
but still more than twice as large as for the (001) sur- 
faces. The unit cell of the 7-plane O-terminated slab has 
the same contents as three bulk unit cells, so the relevant 
surface energy is just 

E,^,f{0) = i (O) - 3ii;buik] , (5) 



E. CaTiOs (Oil) surface charge distributions and 
chemical bondings 

Wc present in Table VIII the calculated MuUiken ef- 
fective charges Q, and their changes AQ with respect to 
the bulk values, for atoms near the surface for the various 
(Oil) surface terminations. 

On the TiO-terminated surface, the charge on the sur- 
face Ti atom is seen to be substantially reduced relative 
to the bulk, while the metal atoms in the third layer lose 
much less charge. The O ions in all layers except the 
central one also have reduced charges, making them less 
negative. The largest charge change (0.232 e) is observed 
for subsurface O atoms, giving a large positive change of 
0.464 e in the charge for that subsurface layer. 

On the Ca-terminated surface, negative changes in the 
charges arc observed for all atoms except for the oxygens 
in the central layer and the Ti atom in the third layer. 
The largest charge changes are for the surface Ca ion 
and the subsurface O ion. The largest overall change in 
a layer charge (— 0.234 e) appears in the subsurface layer 
as well. 

For the O-terminated surface, the negative charge on 
the surface oxygen is very strongly decreased. Corre- 
spondingly, the second layer becomes substantially more 
negative (overall change — 0.177e), with the change com- 
ing mostly on the Ti atom. The total charge density on 
the third layer is almost unchanged. Negative changes in 
charge are observed on all central layer atoms, leading to 
a total charge change of — 0.055 e in that layer. 

The interatomic bond populations for the three termi- 
nations of the (Oil) surface are given in Table IX. The 
major effect observed here is a strong increase of the Ti-0 
chemical bonding near the TiO- and O-terminated sur- 
face as compared to bulk (0.084 e) or to what was found 
on the Ti02-terminatcd (001) surface (0.114e). For the 
O-terminated surface, the 0(I)-Ti(II) bond population is 
about twice as large as in the bulk, and about half again 
as large as at the Ti02-terminated (001) surface. For 
the TiO-terminated (Oil) surface, the Ti-0 bond pop- 
ulations are larger in the direction perpendicular to the 
surface (0.186 e) than in the plane (0.128 e). 
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TABLE VIII: Calculated MuUiken atomic charges Q, and 
their changes AQ with respect to the bulk, in e, for the three 
CaTiOs (Oil) surface terminations. For reference, the bulk 
values are 2.330 e (Ti), -1.371 e (O), and 1.782 e (Ca). 



Atom (layer) Q AQ 

TiO-terminated surface 

Ti(I) 2.204 -0.126 

0(I) -1.290 0.081 

0(II) -1.139 0.232 

Ca(III) 1.733 -0.049 

Ti(III) 2.309 -0.021 

O(III) -1.302 0.069 

0(IV) -1.375 -0.004 

Ca-terminated surface 

Ca(I) 1.676 -0.106 

0(II) -1.488 -0.117 

Ca(III) 1.781 -0.001 

Ti(III) 2.334 0.004 

O(III) -1.452 -0.081 

0(IV) -1.363 0.008 

O-terminated surface 

0(1) -1.139 0.232 

Ca(II) 1.751 -0.031 

Ti(II) 2.235 -0.095 

0(11) -1.422 -0.051 

O(III) -1.359 0.012 

Ca(IV) 1.774 -0.008 

Ti(IV) 2.310 -0.020 

0(IV) -1.398 -0.027 



IV. CONCLUSIONS 

According to the results of our ah initio hybrid B3PW 
calculations, all of the upper-layer atoms for the Ti02- 
and CaO-terminated CaTiOa (001) surfaces relax in- 
wards, while outward relaxations of all atoms in the 
second layer are found at both kinds of (001) termina- 
tions. These results are typical for other technologically 
important ABO3 perovskites such as BaTiOs, PbTiOa, 
and BaZrOsi^i^i^i^ However, they contrast with the 
only previous ah initio study of CaTiOs (001) surfaces 
by Wang et air- where the authors found that the first- 
layer O atoms relax outwards on the CaO-terminated 
(001) surface. For the Ti02-terminated (001) surface, 
our largest relaxation is on the second-layer atoms, not 
on the first-layer ones, this time in agreement with Wang 
et alM- The stronger relaxation of the second-layer atoms 
compared to the first-layer ones was found by us ear- 
lier also for TiOz-terminated PbTiOg and SrTiOa (001) 
surfacesj^i^ Our calculations of the CaO-terminated 
(001) surface shows a very strong inward relaxation of 
8.31% for the top-layer Ca atoms, in very good quan- 
titative agreement with the inward relaxation of 8.80% 



TABLE IX: The A-B bond populations P (in e) and the rel- 
evant interatomic distances R (in A) for three different (Oil) 
terminations of the CaTiOa surface. Symbols I-IV denote the 
number of each plane enumerated from the surface. The near- 
est neighbor Ti-0 distance in the unrelaxed bulk is 1.926 A. 



Atom A 


Atom B 


P 


R 


TiO-terminated surface 






Ti(I) 


0(1) 


0.128 


1.979 




0(11) 


0.186 


1.752 


0(II) 


Ti(III) 


0.110 


1.935 




Ca(III) 


0.018 


2.769 




O(III) 


-0.024 


2.790 


Ti(III) 


Ca(III) 


0.000 


3.336 




O(III) 


0.100 


1.929 




0(IV) 


0.076 


1.904 


Ca(III) 


O(III) 


0.008 


2.723 




0(IV) 


0.004 


2.672 


O(III) 


0(IV) 


-0.032 


2.653 


Ca-terminated surface 






Ca(I) 


0(11) 


0.006 


2.458 


0(II) 


Ca(III) 


0.012 


2.768 




Ti(III) 


0.072 


1.973 




O(III) 


-0.036 


2.784 


Ca(III) 


O(III) 


0.002 


2.723 




0(IV) 


0.006 


2.705 


Ti(III) 


O(III) 


0.060 


1.926 




Ca(III) 


0.000 


3.335 




0(IV) 


0.084 


1.915 


O(III) 


0(IV) 


-0.064 


2.691 


O-terminated surface 






0(I) 


Ca(II) 


0.028 


2.613 




Ti(II) 


0.162 


1.699 




0(11) 


-0.016 


2.788 


Ca(II) 


0(II) 


-0.006 


2.412 




Ti(II) 


0.002 


3.198 


Ti(II) 


0(II) 


0.086 


1.992 




O(III) 


0.100 


1.764 


0(11) 


O(III) 


0.010 


2.925 


Ca(II) 


0(111) 


0.006 


2.737 


O(III) 


0(IV) 


-0.038 


2.750 




Ti(IV) 


0.062 


1.963 




Ca(IV) 


0.004 


2.677 



found by Wang et alA This inward relaxation of the sur- 
face Ca atoms on the CaO-terminated (001) surface is 
much stronger than was obtained for the AO-terminated 
(001) surfaces of other ABO3 perovskites (A — Sr, Ba, 
Pb, and Zr)i^^>Si>62 

Our calculated surface rumphng of 7.89% for the 
CaO-terminated (001) surface is almost five times larger 
than that of the corresponding Ti02-terminated sur- 
face, and is comparable with the surface rumpling of 
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9.54% obtained for the CaO-terminated surface by Wang 
et alP This rumphng is larger than the rumphngs ob- 
tained in previous ab initio calculations for the AO- 
terminated (001) surfaces of SrTiOs, BaTiOs, BaZrOs, 
and PbTiOs.^'^^^'^'^^ 

Our calculations predict a compression of the inter- 
layer distance between first and second planes, and an 
expansion between second and third planes, for the (GDI) 
surfaces. Our value for Ac?i2 of —9.43% on the CaO- 
terminated (001) surface is in a reasonable agreement 
with the result of —11.43% obtained by Wang et a.^ and 
is larger than the corresponding value for 3, BaTiOa, 
BaZrOs, and PbTiOg (001) surfaces^i^i2ii32^^^ 
for experimental confirmation of these results, we are 
unfortunately unaware of experimental measurements of 
Adi2 and Adas for the CaTiOs (001) surfaces. More- 
over, for the case of the SrO-terminated SrTiOa (001) 
surface, existing LEED^^ and RHEED^^ experiments ac- 
tually contradict each other regarding the sign of A(ii2. 
In view of the absence of clear experimental determina- 
tions of these parameters, therefore, the first-principles 
calculations are a particularly important tool for under- 
standing the surface properties. 

Turning now to the CaTiOa (Oil) surfaces, we found 
that the inward relaxation of the upper-layer metal 
atom on the TiO-terminated (Oil) surface (Ti displace- 
ment of 7.14%) is smaller than on the CaO-terminated 
(001) surface (Ca displacement of 8.31%), in contrast to 
what was found for the SrTiOg, BaTiOg, PbTiOg, and 
BaZrOa surfaces However, the inward re- 
laxation by 16.05% of the upper-layer Ca atom on the 
Ca-terminated (Oil) surface is about twice as large as 
the inward relaxations of surface atoms obtained on the 
CaO-terminated (001) surface. Our calculated atomic 
displacements in the third plane from the surface for 
the Ca, TiO, and 0-terminated (Oil) surfaces are still 
substantial. Our calculated surface rumpling s for the 
TiO-terminated (Oil) surface is approximately 1.5 times 
larger than that of the CaO-terminated (001) surface, 
and many times times larger than that of the Ti02- 
terminated (001) surface. Also, our ah initio calculations 
predict a compression of the interlayer distance A(ii2 and 
an expansion of Ac?23 for the TiO- and 0-terminated 
(Oil) surfaces. This behavior seems to be obeyed by all 
previous calculations of relaxations at (001) ABO3 per- 
ovskite surfacesii2£i2ii22iSiiS^; we can conclude that this 
effect may be a general rule, requiring further experimen- 
tal studies and confirmation. 

A comparison of our ah initio B3PW calculations on 
the TiO-terminated CaTiOa (Oil) surface with the pre- 
vious ah initio calculation performed by Zhang et ali^ 
shows that the atomic displacement directions almost al- 
ways coincide, the only exception being the small third- 
layer Ti-atom inward relaxation of —0.78% in our cal- 
culation compared with an outward one of 0.28% in 
theirs. The displacement magnitudes are generally com- 
parable in the two studies, leading to an excellent agree- 
ment for the TiO-terminated (Oil) surface rumplings 



(11.81% in our calculations vs. 12.10% in theirs). For 
the Ca-terminated (Oil) surface, our inward relaxation 
magnitude of 16.05% for the upper-layer Ca atom, the 
largest of all atoms on all of the studied (Oil) termina- 
tions, is in excellent agreement with the value of 15.37% 
obtained in Ref. 6]. This largest displacement of the 
surface A atom on the A-terminated (Oil) surface is 
was also obtained for the SrTiOs, BaTiOa, PbTiOs, 
and BaZrOg cases^i^i^i 132^2^162 j^g^. ^-^ey did for 

the TiO-terminated (Oil) surface, our relaxation direc- 
tions for the Ca-terminated surface almost all coincide 
with those obtained previously)^ the only exception be- 
ing again the displacement direction of the third-layer 
Ti atom. We find that this atom moves slightly inwards 
by 0.37%, whereas the previous work obtain an outward 
relaxation of 0.89%.^ 

For the 0-terminated (Oil) surface, in most cases 
our calculated displacement directions are in qualitative 
agreement with the results of Ref. @. In some cases, 
as for example for the second layer Ti and O atom dis- 
placements in the direction along the surface, our cal- 
culated displacement magnitudes for Ti (4.70%) and for 
O (8.05%) are in an excellent agreement with the cor- 
responding results (4.53% and 8.06% respectively) of 
Zhang et al.^ However, in many cases, our calculated 
displacement magnitude is smaller than that calculated 
in Ref. Most disturbingly, in three cases there are 
also some qualitative differences between our results and 
those of Zhang et al^ Specifically, the second-layer Ca 
and Ti atoms move substantially inwards in our calcu- 
lations, but outwards in Ref. 0], and the third- layer O 
atoms move in opposite directions in the two calculations. 

As for the surface energies, we find that both the CaO- 
and Ti02-terminated (001) surfaces are about equally 
favorable, with surface energies of 0.94 and 1.13 eV re- 
spectively. These values are in excellent agreement with 
the corresponding values of 0.824 and 1.021 eV respec- 
tively as computed by Zhang et al. in Ref. G*]. In con- 
trast, we see very large differences in surface energies 
on the (Oil) surfaces. Our lowest-energy (Oil) sur- 
face is the 0-terminated one at 1.86 eV, with the Ca- 
terminated surface just behind at 1.91 eV, and the TiO- 
terminated surface is very unfavorable at 3.13 eV. These 
are all much larger, by about a factor of two or more, 
than for the (001) surfaces. This is the same ordering 
of (Oil) surface energies as was obtained by Zhang et 
al., but these authors obtained quite different values of 
0.837, 1.671, and 2.180eV for the 0-, Ca-, and TiO- 
terminated (Oil) surfaces, respectivelyjS The values for 
the Ca- and TiO-terminated surface energies are only 
modestly smaller than ours, but the value for the O- 
terminated (Oil) surface energy presents a clear disagree- 
ment with the present work, being more than twice as 
small as ours. In fact, according to their work, the O- 
terminated (Oil) surface is even lower in energy than 
the Ti02-terminated (001) surface, and about equal to 
that of the CaO-terminated (001) surface. In this re- 
spect, their result contrasts not only with our result for 
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CaTiOs, but with all previous ab initio and shell-model 
calculations dealing with SrTiOa, BaTiOa, PbTiOa, and 
BaZrOg (001) and (Oil) surface energies ,-^'2ii23^i32 42 
where the (001) surface energies are always smaller than 
the (Oil) surface energies. 

We do not understand the reason for this discrepancy. 
We have carried out test calculations of the cleavage en- 
ergies of the three (Oil) surfaces using the PBE-GGA 
exchange-correlation functional^"* used by Zhang et al., 
but within the CRYSTAL-2003 code package, and we 
find cleavage energies that are only about 15-25% larger 
than theirs. The drastic difference, then, must be in the 
relaxation energy of the Ca-terminated surface, which 
is — 1.83cV in our calculation and — 2.70eV in theirs. 
We also did a test calculation of the energy of a Ca- 
terminated (Oil) slab in which the surface atoms were 
placed by hand at the coordinates reported for this sur- 
face in Ref. ^6] , and found that the energy was even higher 
than the energy of the unrelaxed structure. Clearly these 
discrepancies call for further exploration. 

Our ab initio calculations indicate a considerable in- 
crease in the Ti-0 bond covalency near the TiO- and O- 



terminated (Oil) surfaces, as well as the Ti02-terminated 
(001) surface. The Ti-0 bond covalency at the TiO- 
terminated (Oil) surface (0.128e) is much larger than 
that for the Ti02-terminated (001) surface (0.114e) or 
in bulk CaTiOs (0.084 e). The Ti-0 bond populations 
on the TiO-terminated (Oil) surface are much larger in 
the direction perpendicular to the surface than in the 
plane (0.186 vs. 0.128 e). Our calculated increase of the 
Ti-O bond covalency near the (Oil) surface, is in agree- 
ment with the resonant photoemission experiments.— 
This should have an impact on the electronic struc- 
ture of surface defects (e.g., F centers), as well as on 
the adsorption and surface diffusion of atoms and small 
molecules relevant for catalysis. 
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